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Abstract. The evolution of an initially stellar dipole type magnetosphere interact¬ 
ing with an accretion disk is investigated using numerical ideal MHD simulations. 
The simula tions follow several 1000 Kepler ian periods of the inner disk (for animated 
movies see ^ttp ://www.aip.de/^cfendt ). Our model prescribes a Keplerian disk 
around a rotating star as a fixed boundary condition. The initial magnetic field 
distribution remains frozen into the star and the disk. The mass flow rate into the 
corona is fixed for both components. 

The initial dipole type magnetic field develops into a spherically radial outflow 
pattern with two main components — a disk wind and a stellar wind — both evolving 
into a quasi-stationary hnal state. A neutral held line divides both components, along 
which small plasmoids are ejected in irregular time intervals. The half opening angle 
of the stellar wind cone varies from 30° to 55° depending on the ratio of the mass 
flow rates of disk wind and stellar wind. The maximum speed of the outhow is about 
the Keplerian speed at the inner disk radius. 

An axial jet forms during the hrst decades of rotations. However, this feature 
does not survive on the very long time scale and a pressure driven low velocity how 
along the axis evolves. Within a cone of 15° along the axis the formation of knots 
may be observed if the stellar wind is weak. With the chosen mass how rates and 
held strength we see almost no indication for a how self-collimation. This is due to 
the weak net poloidal electric current in the magnetosphere which is in difference to 
typical jet models. 


1. Introduction 


A stellar dipole magnetic field surrounded by an accretion disk is the 
model scenario for a variety of astrophysical sources — classical T Tauri 
stars, cataclysmic variables or high mass X-ray binaries. Some of them 
exhibit Doppler shifted emission lines indicating wind motion. Highly 
collimated jets are observed from young stellar objects and X-ray bina¬ 
ries. In general, magnetic fields are thought to play the leading role for 
jet acceleration and collimation (Blandford &: Payne 1982; Camenzind 
1990; Shu et al. 1994; Fendt et al. 1995). 

Several papers numerically consider the evolution of a stellar mag¬ 
netic dipole interacting with a diffusive accretion disk (Hayashi et al. 
1996; Miller &: Stone 1997; Goodson et al. 1997, 1999). In these papers 
a collapse of the inner disk is indicated. The inward accretion flow 
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develops a shock near the star. The stream becomes deflected resulting 
in a high-speed flow in axial direction. The results of Goodson et al. 
are especially interesting as combining a huge spatial scale with high 
spatial resolution near the star. However, to our understanding, it is not 
clear, how the initial condition (a standard a-viscosity disk) is actually 
evolving in their code without physical viscosity. 

We emphasize that time-dependent simulations lasting only a short 
time period strongly depend on the initial condition. The simulation 
of jet magnetosphere over many rotational periods is essential. The 
observed kinematic time scale of protostellar jets can be as large as 
lO^yrs (or 5 x 10^ stellar rotations). Thus, we follow another approach 
for the simulation of magnetized winds from accretion disks and con¬ 
sider the accretion disk “only” as a boundary condition. Since the disk 
structure itself is not treated, such simulations may last over hundreds 
of Keplerian periods. For a certain initial magnetic field distribution, 
a stationary state self-collimating jet flow can be obtained (Ouyed & 
Pudritz 1997, Fendt & Cemeljic 2002). 

In our project (Fendt & Elstner 1999, 2000), we are essentially in¬ 
terested in the evolution of the ideal magnetohydro dynamic (MHD) 
magnetosphere and the formation of winds and jets and not in the 
evolution of the disk structure itself. Therefore, we do not include 
magnetic diffusivity into our simulations. We follow the approach of 
Ouyed & Pudritz (1997) in combination with dipole magnetic held 
initial condition. Using the ZEUS-3D code (Stone & Norman 1992a,b) 
in the axisymmetry option we solve the time-dependent ideal MHD 
equations, 

^ + y ■ (pv) = 0, ^-Vx(FxH) = 0, V-.B = 0, 


+ V{P + Pa) + -j X B = 0, 

held, V the velocity, p the gas density, P 
the gas pressure, j = V x BjAir the electric current density, and 
the gravitational potential. We apply a polytropic equation of state 
(polytropic index 7 = 5/3) and do not solve the energy equation. 
Instead, the internal energy is defined given e = p /(7 — 1). Similar 
to Ouyed & Pudritz (1997), we introduce an Alfven wave turbulent 
magnetic pressure. Pa = P/Pt, with constant Pt- 

The main parameters of our simulation are the plasma beta just 
above the inner disk radius r*. Pi = Si^Pi/B^, and the Mach number of 
the rotating gas, 5i = Piv\^/Pi, where VK,i = a/ GM /. 


dv 

m 


-I- (F- V) F 


where B is the magnetic 
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2. The model numerical realization 

Our model setup represents a central star surrounded by a Keplerian 
disk with a gap in between. 

The initial coronal density distribution is in hydrostatic equilibrium. 
The gravitational potential is not smoothed. We choose the initial 
field distribution of a force-free, current-free stellar magnetic dipole, 
deformed by the effect of “dragging” in the disk. Therefore the poloidal 
field is inclined towards the disk surface. Since the boundary condition 
for the poloidal magnetic field along the inflow boundary is fixed, the 
magnetic flux from the star and disk is conserved. The disk toroidal 
magnetic held is also force-free, = //j/r for r >r\ with jii = 

Hydrodynamic boundary conditions are ‘inhow’ along the r-axis, ‘re- 
hecting’ along the symmetry axis and ‘outhow’ along the outer bound¬ 
aries. The inhow parameters into the corona are dehned with respect to 
the three different boundary regions - star, gap and disk. The stellar 
wind boundary condition is motivated by the fact that stellar winds 
are indeed observed. The ratio of mass how rates in the two outhow 
components also governs the structure of the how. 


3. Results and discussion 

In the following we discuss the results of our simulations. For details see 
Fendt & Elstner (1999, 2000). As a general behavior, the initial dipole 
type structure of the magnetic held disappears on spatial scales larger 
than the inner disk radius and a two component wind structure - a 
disk wind and a stellar wind - evolves. On the very long time scale we 
hnd a quasi-stationary hnal state of a spherically radial mass outhow. 

The stellar rotational period is chosen as H* = {vK,i/f’i) = 1. Thus, 
the magnetospheric co-rotation radius is at r*. Other parameters are 
f5i = 1.0, 5i = 100, ^ = —1.0, and a stellar radius r* = 0.5. A numerical 
grid of 250^ elements is used for a box of physical size 20 x 20rj. The 
stellar wind mass how rate is rather large, M*/Md = 2. 

3.1. The first evolutionary stages 

During the hrst evolutionary stages the magnetospheric structure is 
characterized by the following features. The winding-up of the dipolar 
poloidal held, the formation of a neutral field line, a transient axial jet 
feature, a two component outhow consisting of a stellar wind and a disk 
wind. 

The winding-up process of poloidal magnetic held due to differential 
rotation between star and disk and the static initial corona induces a 
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Figure 1. Long-term evolution of a axisymmetric dipolar magnetosphere. Shown is 
density (grey scale) and poloidal magnetic field (white lines) for 0, 10, 25, 50, 100, 
250, 500, 950, and 2700 rotations. The symmetry axis points upwards. 


toroidal field with a positive (negative) sign along the field lines outside 
(inside) the slowly emerging neutral field line. 

During the first 50 rotations a jet feature evolves along the rotational 
axis with a pattern velocity of about 0.2 UK,i- Such an axial jet is known 
as a characteristic result of MHD simulations performed in the recent 
literature (Hayashi et ah 1997; Goodson et al. 1997, 1999; Kudoh et 
al. 1998). However, we find that the formation of this feature results 
from the adjustment of the initially hydrostatic state to a new dynamic 
equilibrium and disappears on the long time scale. 

The disk wind accelerates rapidly from the low injection speed to 
fractions of the Keplerian speed. The flow starts already super Alfvenic 
due to the weak dipolar field in the disk. Thus, magneto-centrifugal 
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Figure 2. Long-term evolution of a dipolar magnetosphere. The velocity vectors for 
the time step of 30 and 2500 rotations clearly show the two flow components and 
the weakening of the axial jet. The components are separated by the neutral field 
line. 


acceleration along the inclined dipole type field lines of the initial 
magnetic held is not the acceleration mechanism. The acceleration 
mechanism is mainly due to the centrifugal force on the disk mat¬ 
ter reaching the non-rotating corona. Higher above the disk also the 
Lorentz force contributes to the acceleration. 

The rotating stellar magnetosphere generates a stellar wind. Due 
to the strong stellar held, the how starts sub-Alfvenic. It is initially 
magneto-centrifugally driven with a roughly spherical Alfven surface. 

3.2. The long-term evolution 

The total mass how rate into the corona determines how fast the how 
will establish a (quasi-)stationary state. First, the outhow evolution 
is highly time-variable and turbulent. After relaxation of the MHD 
conhguration from the initial magnetohydrostatic state into a new 
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dynamical equilibrium, we finally observe a two-component outflow 
from disk and star distributed smoothly over the whole hemisphere 
and moving in spherically radial direction. After about 2000 rotations 
a quasi-stationary outflow is established over the whole grid (Fig. 1). 
However, small scale instabilities can be observed along the neutral field 
line separating stellar and disk wind and along the symmetry axis. 

Figure 2 shows the poloidal velocity vectors. High velocities (> UK,i) 
are only observed far from the axis. The asymptotic speed is about 
1.5 UK,i for both components. The initial axial jet feature disappears. 
The axial flow moves with 0.2uK,i- 

The axial blobs (or rather tori) generated in this simulation run 
move with pattern speed of about 0.1 the Keplerian speed at r*. We 
emphasize that the knot size and time scale of knot formation in our 
simulation is far from the jet knots observed in protostellar jets 

The quasi-stationary two-component outflow obtained in our simu¬ 
lations show almost no indication for collimation. This is in agreement 
with the analysis of Heyvaerts & Norman (1989) who have shown that 
only jets carrying a net poloidal current will collimate to a cylindrical 
shape. In our case we have an initially dipolar magnetosphere and 
the final state of a spherically radial outflow enclosing a neutral line 
with a poloidal magnetic field reversal. The toroidal field reversal im¬ 
plies a reversal also of the electric current and, thus, only a weak net 
poloidal current. Thus, the flow self-collimation naturally obtained in a 
monotonous magnetic flux distribution (Ouyed & Pudritz 1997, Fendt 
& Cemeljic 2002) cannot be achieved. 
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